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Motivated by the possibility that nonets of scalar mesons might be described as mixtures of "two 
quark" and "four quark" components, we further study a toy model in which corresponding chiral 
nonets (containing also the pseudoscalar partners) interact with each other. Although the "two 
quark" and "four quark" chiral fields transform identically under SU(3)lX SU(3)a transformations 
they transform differently under the U(1)a transformation which essentially counts total (quark + 
antiquark) content of the mesons. To implement this we formulate an effective Lagrangian which 
mocks up the U(1)a behavior of the underlying QCD. We derive generating equations which yield 
l/") ' Ward identity type relations based only on the assumed symmetry structure. This is applied to the 

mass spectrum of the low lying pseudoscalars and scalars. as well as their "excitations". Assuming 
isotopic spin invariance, it is possible to disentangle the amount of "two quark" vs. "four quark" 
content in the pseudoscalar it, K, n type states and in the scalar k type states. It is found that a small 
"four quark" content in the lightest pseudoscalars is consistent with a large "four quark" content 
in the lightest of the scalar k mesons. The present toy model also allows one to easily estimate the 
strength of a "four quark" vacuum condensate. There seems to be a rich and interesting structure. 
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I. INTRODUCTION 



> 
O 

The last few years have seen a renewal of interest Q- in the low energy scalar sector of QCD. Many physicists 
^sO \ now believe in the existence of the light, broad I = J = resonance, sigma in the 500-600 MeV region as well as a 
light broad I,J= 1/2, resonance, kappa in the 700-900 MeV region. Together with the well established /o(980) and 
ao(980) scalar resonances, these comprise a putative nonet of "elementary particles". Furthermore, this nonet seems 
likely to have a quark structure like qqqq rather than the conventional qq [3l| . This of course raises the question of 
,_SZh ' where are the conventional qq p-wave scalars expected in the quark model. Arguments have been given 32] that the 
experimental data are better fit when the two scalar nonets mix with each other and the resulting "level repulsion" , 
O,' pushes the conventional scalars to higher masses than otherwise expected. 

In order to further explore the feature of mixing between qq type and qqqq type states it seems interesting to 
1 consider a linear SU(3) x SU(3) sigma model which contains also the pseudoscalar nonet partners of these two scalar 
nonets. Parenthetically, we remark that while the non- linear sigma model [33l l34| and its extension to the chiral 
perturbation theory program [35| are often more efficient for systematic calculations, linear sigma models have a very 
long history of furnishing important insights into the nature of strong hadron dynamics. The SU(2) linear sigma 
model was first given in ref. I33|. It was used as a basis for understanding the current algebra treatment of hit 
scattering near threshold in ref. |36j| . The SU(3) version was given in the first of ref. [37| ■ A detailed application to 
the low energy pseudoscalar mass spectrum was given [38[ before QCD in which, among other things, it was shown 
how a U(1)a violating term natural in the SU(3) model could solve the rf problem. Such a term was later discovered 
to arise from instanton effects [3|| . The connection was pointed out in ref. [ifj and emphasized by 't Hooftlllj • 

The model containing two different chiral nonets to be discussed here was proposed in section V of ref. |2!j and an 
initial treatment, neglecting flavor symmetry breaking, was given. A discussion, taking the flavor symmetry breaking 
into account has very recently been presented in ref. [42|. Actually, it turns out that the model is very complicated 
since many different terms can be included and various assumptions about the nature of the symmetry breaking can 
be made. In this paper we will set up the formalism for treating consequences of the model which hold (at tree level) 
just due to the symmetry structure of the model and will give a numerical treatment using what might be the simplest 
choice of symmetry breaking terms. 

Section II begins with a review of the flavor transformation properties of the two chiral nonet fields, M and M' 
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which are used in the model. Each contains nine pseudoscalar and nine scalar fields. Under chiral SU(3)lX SU(3)a 
transformations both fields transform in an identical manner. Thus a chiral Lagrangian which respects only this 
symmetry cannot directly distinguish between a "two quark" (i.e. qq) or a "four quark" scalar, for example. However, 
it is noted again that the U(1) J 4 transformation actually counts the number of quarks in these mesons and provides a 
way to distinguish them. In order to make use of this, the Lagrangian should of course be set up appropriately. We 
implement this by requiring that the Lagrangian mock up the anomalous U(1)a equation of the underlying QCD and 
that the analogs of the quark mass terms also mock up the U(1)a transformation properties of the quark mass terms 
in the underlying theory. A reasonable initial thought on which terms to include in the Lagrangian is to restrict it to 
be renormalizable. It is noted, with details in Appendix A, that the renormalizable M — M' Lagrangian has however 
very many more terms than does the renormalizable single M Lagrangian. An alternate way, which still satisfies 
generality, is to consider any number of terms, renormalizable or not, and just use the information which follows from 
the symmetry behavior of the Lagrangian. 

In order to exploit this symmetry information we derive, in section III, vector type and axial vector type "generating 
equations" for the model. These can be differentiated with respect to the fields to yield many tree level Ward identities 
which arc independent of the number of symmetric terms included in the Lagrangian. In addition to the analog of 
"two quark" condensates which occur in the single M model, the present model also brings "four quark" condensates 
into the picture. 

In section IV, we derive predictions for the mass spectrum which follow from this symmetry approach. The 
characteristic feature is mixing between "two quark" and "four quark" mesons with the same quantum numbers. 
Assuming isospin invariance, predictions are made for the tt — it' mixing sector, the K — K' mixing sector, the strange 
scalar k — k' mixing sector and the sector involving mixing of the four isocalar pseudoscalars (77 type particles) . It is 
shown how to formulate the first three of these mixing sectors in a parallel and economical way. 

In section V, the mass spectrum relations are compared with experiment. First the three 2x2 mixing sectors are 
treated. The inputs are taken to be the six masses of the well known and not so well known particles, the pion and 
kaon decay constants and a model parameter denoted x n , which is the squared mass of the unmixed (or "bare" ) pion. 
These are enough to determine all the relevant parameters of these three systems. The pseudoscalar mixing is very 
sensitively dependent on as it increases from the experimental value, the four quark components of the pion 
and the kaon increase. On the other hand, the scalar k has a large four quark component. This feature thus provides 
some support for a more exotic structure of the low lying scalars. Another interesting feature of the present model, 
discussed in this section, is that it permits one to estimate the strength of a four quark vacuum condensate. Finally, 
section V contains a brief summary, the connection with other results on the same model and directions for future 
work. 



II. SYMMETRIES AND LAGRANGIAN 



First, let us briefly review 29] the fields of the model and their transformation properties. The schematic structure 
for the matrix M(x) realizing a qq composite in terms of quark fields q a A{x) can be written 

M b a = ((?bA )t 74 i±^ 9QA , (1) 

where a and A are respectively flavor and color indices. Our convention for matrix notation is M h a — > M a f,. Then M 
transforms under chiral SU(3)lX SU(3)r as 

M -» U L MU R , (2) 

where Ul and Ur are unitary, unimodular matrices associated with the transformations on the left handed 
(ql = 5 (1 + 75) q) and right handed (qn — | (1 — 75) q) quark projections. For the discrete transformations charge 
congugation C and parity P one verifies 

C : M — > M T , P : M(x) -► A/ f (-x). (3) 

The \J(1)a transformation acts as q a L —> e ll, q a L, q a R — * e~ lv qaR and results in: 

M -> e 2iv M. (4) 

Next, consider the qqqq type fields. One interesting model |43j postulates that the light scalars are "molecules" made 
out of two pseudoscalar mesons. The chiral realization of this picture would result in the following schematic structure: 

M^ b = e acd e^ \M%{M% (5) 
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One can verify that transforms exactly in the same way as M under SU(3)lX SU(3)#, C and P. Under XJ(1)a 
it transforms as 

M (2) ^ e -4fe M (2) 5 ( 6 ) 

which differs from Eq. Q. Another interesting approach [3l| to explaining the light scalar mesons was formulated 
by Jaffe in the framework of the MIT bag model. It was observed that the spin-spin (hyperfine) piece of the one 
gluon exchange interaction between quarks gives an exceptionally strong binding to an s-wave qqqq scalar state. The 
scalar states of this type may be formally written as bound states of a "dual quark" and "dual antiquark" . There are 
two possibilities if the dual antiquark is required to belong to a 3 representation of flavor SU(3). In the first case it 
belongs to a 3 of color and is a spin singlet. This has the schematic chiral realization, 

rgE gab EAB T r<-X ^ 5 - 

L =e € q a A^ — 2 — qhBl 

m E — ,ffab BAB T n -l 1 ~ 75 /_s 

W - e y e q aA C — - — q bB , (7) 

where C is the charge conjugation matrix of the Dirac theory. A suitable form for the M matrix is: 

AfW = {U A )^R^ A . (8) 

can be seen to transform in the same way as under SU(3)lX SU(3)a, C, P and U(1)a- In the second case 
the dual antiquark belongs to a 6 representation of color and has spin 1. It has the corresponding schematic chiral 
realization: 



t9 — t9 - 9 a0 T _ j" 75 

nv.AB — ^u,v,BA — e 1aA^ °V o qbB ' 



RLab = Ku.ba = e gab q T aA C- l a^^-^q bB , (9) 



J liv,AB ~ ^nv,BA ~ c VaA^ u y.v 2 

D9 _ n9 

l fiv,AB ~ "-iiv,BA 

where — A [7^,7^]- This choice leads to an M matrix 

^'^W^- (10) 

where the dagger operation includes a factor (_i)' 5 f 4 +' 5 >' 4 _ j\f (4) a } so transforms like M( 2 ) and M^ 3 ) under all of 
SU(3)lX SU(3)i{, C, P and U(1)a- The specific form favored by the MIT bag model calculation actually corresponds 
to a particular linear combination of and M^ 4 \ Furthermore one can verify that in Eq. (JSJ is related by a 
Fierz transformation to a linear combination of and M^ 4 \ Thus only two of M^ 2 \ and are linearly 

independent. In any event, at the present effective Lagrangian level, there are no quantum numbers to distinguish 
M<?\ Af( 3 ), and from each other so we may as well just denote an arbitrary linear combination of them to be 

our qqqq field, M' . Note that M and M' are distinguished from each other by their different U(1)a transformation 
properties. These fields may be decomposed into hermitian scalar (S) and pseudoscalar (<j>) nonets as, 

M = S + i<j>, 

M' = S' + ic/)'. (11) 

We will be interested in the situation where non-zero vacuum values of the diagonal components of S and S' may 
exist. These will be denoted by, 

(S b a )=a a 5 b a , (S' a b ) = (3 a S b a . (12) 

In the iso-spin invariant limit, ot\ — and /3i — /3a while in the SU(3) invariant limit, ot\ — ct2 — ol^ and fli = 02 = Pz- 
The Lagrangian density which defines our model is 

C = -~Tr (d^Md^) - ^Tr (d^M'd^) - V (M, M 1 ) - V SB , (13) 

where Vq(M, M') stands for a general function made from SU(3)z,x SU(3)# (but not necessarily U(l)yt) invariants 
formed out of M and M' . Furthermore V$b is taken to be a flavor symmetry breaking term which should mock up 
the quark mass terms which perform this function in the fundamental QCD Lagrangian. Other physical particles 
(including glueballs) could be added for more realism, but Eq. Q13f) is already quite complicated. 
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To get an initial indication of what is happening in this kind of model the drastically simplified case where the 
quark mass effective term, V$b is absent and where Vq is simply given by: 

Vq = -c 2 Tr (M M f ) + c 4 Tr (MM^MM ] ) + d 2 Tr (M'M*) + eTr (MM n + M'M ] ) , (f 4) 

was treated in sec.V of ref. [2j|. Here C2, C4 and d 2 are positive real constants. The M matrix field is chosen to 
have a wrong sign mass term so that there will be spontaneous breakdown of chiral symmetry. A pseudoscalar octet 
is thus massless. The mixing between the M and M' is controlled by the parameter e. The first feature found for 
this simplified model was that the analog, (S'") of the qqqq condensate in QCD acquired a small non-zero value 
due to the mixing between S and S' . The main question is the level ordering. Since the light pseudoscalars (e.g. 
7r + = <f>\) arc naturally identified, before mixing, with the qq field M, one wonders whether the two quark rather 
than the four quark scalars aren't the lightest ones. It was found however that it is natural (but not unique) in the 
model to have the energy level pattern in ascending order- pseudoscalar Nambu-Goldstone boson with primarily qq 
structure, scalar with primarily qqqq structure, pseudoscalar with primarily qqqq structure and scalar with primarily 
qq structure. These refer to degenerate octets which are each mixtures of M and M' states. This seems to be similar 
to the expected experimental pattern and gives us some motivation to proceed further. 

The next question is what terms to include in the Lagrangian Eq. (|f 31) . A natural first attempt would be to 
consider a renormalizable model in which Vq contains all the SU(3) x SU(3) invariant terms up to four powers of the 
fields. These are listed in Appendix A. It is seen that there are 21 terms of this type. This is a rather large number 
and while not impossible to handle suggests trying another tack. We will just allow Vb to contain all possible terms 
which are SU(3) l x SU(3) # symmetric and use the information provided by this symmetry. This is more general and 
also allows for non-renormalizable terms. The price to be paid is that we only get information which follows just from 
the symmetry structure. In an earlier treatment |38| of the single chiral nonet case, it was found that the results 
obtained were essentially those which could be obtained from the "current algebra" approach. Furthermore, we will 
try to make use of the fact that M and M' have different U(1)a transformation properties. We thus demand that the 
Lagrangian without Vsb mock up the anomalous U(1)a equation of QCD, 

SC = G, (15) 

where 8 denotes the axial U(l) variation and G is proportional to the product of the QCD field strength tensor and 
its dual. This can be achieved by making all of the terms in Vb, except for a limited number, XJ(1)a invariant. The 
special terms will be constructed to satisfy Eq. 11511 . An example of a term which is not U(1)a invariant is the mixing 
term used in the simplified model above: Tr (M'M^) + h.c. However a mixing term of the type: 

e abc e def M%M b e M' f c + h.c. (16) 

is U(l) J 4 invariant and hence possibly the most important one. 

An SU(3)lx SU(3)h invariant but not U(l)^ invariant term which mocks up Eq. i(T5|) can be seen 01 to be 

iG detM 
£ ~ = ^ ln( detMt } - 

Here, G is being formally considered as an effective pseudoscalar glueball field in the effective Lagrangian. To get an 
77' (958) mass term in the effective lagrangian framework one can [44| include a wrong sign mass term for G: cG 2 /2 
in the Lagrangian which of course does not change the flavor symmetry structure. Then integrating out G yields the 
effective 7/(960) mass term: 

. detM , , 

where C3 = — l/(288c). The nature of this term becomes more apparent when one goes to the non- linear realization 
where M — > ajexp (i(j>/a.\). For the present paper we shall consider this to be the only SU(3)iX SU(3)/j invariant 
but not \J(\)a invariant term. However, it is not at all unique when we consider a model with two chiral nonets. For 
example one can also include something like the non- U(l)^ invariant mixing term Tr [M 1 M') + h.c. by writing a 
candidate Lagrangian piece: 

iG, , .detM, , .Tr (MM*). , . 

12 [7lln Wt } + 72hl ( Tr(MW' (M) 

and proceeding as above. In order to properly mock up the anomaly in this case it is necessary (4^ that the real 
numbers 71 and 72 satisfy 

71 +72 = 1- (20) 



5 



The generalization to more than two such terms is evident. It may be noted that the M — M' mixing term resulting 
from Eq. I|19l) mixes only the pseudoscalar fields and not the scalar ones. 

Finally, let us consider the flavor symmetry breaking terms. To get more restrictions, we assume that such a term 
should mock up both the SU(3)^x SU(3)a and U(l)^ transformation properties of the quark mass terms in the 
fundamental QCD Lagrangian. It is convenient to introduce a diagonal matrix, 

A = diag(A 1 ,A 2 ,A 3 ), (21) 

which is proportional to the diagonal matrix made from the three light quark masses, diag{m Ul rrid, m s ) (See |4f)j | for 
further details). Then, from Eq. (Q, we note an obvious choice for a flavor symmetry breaking term, 

V SB = -Tr [A(M + Aft)] = -2Tr (AS), (22) 

which transforms like (3,3*) + (3*, 3) under SU(3)z,x SU(3)a. Under the U(1)a transformation of Eq. J3J, it goes 
to — e 2w Ti (AM) +h.c. Note that the similar simple possibility, — 2Tr (AS") does not correctly mock up the U(l)^ 
transformation property of the QCD mass term. However Eq. I)22|l is not at all unique in correctly mocking up the 
quark mass term. An interesting term which does mock up the quark mass term also involves mixing and has the 
form, 

e abc e de fA a d M b e M' c f + h.c. (23) 

This term mixes both scalars and pseudoscalars but with opposite signs. 

For what follows, it is convenient to record the behaviors of the fields under infinitesimal transformations. Let us 
write the infinitesimal vector (L+R) and axial vector (L-R) transformations of cf> and S as, 

5 v <f> = [E v ,(f>], S A (f) = —i[E A ,S]+, 
6 V S = [E V ,S], S A S = i[E ' A , 4>] + . (24) 

Here, unitarity demands that the infinitesimal matrices obey, 

El = -E v , E\ = -E A . (25) 

If we demand that the transformations be unimodular, so that the U(l)^ transformation is not included (the U(l)y 
transformation is trivial for mesons), we should also impose Tr (E A ) = 0. However we will not do this so the effects 
of U(1)a will also be included. The transformation properties of the qqqq type fields are: 

Sytf = [E v , <j>'], 6 A <t>' = —i[E A , S'}+ + 2z5'Tr (E A ), 

S V S' = [E v , S'}, S A S' = i[E A , (/>']+ - 2z0'Tr(£ A ). (26) 

The extra terms for the axial transformations reflect the different U(1)a transformation properties of M and M'. 



III. GENERATING EQUATIONS 

We shall consider, in this paper, tree level predictions for the Lagrangian of Ea. (|13|l in which the only U(1)a 
violating term in Vq is that of Eci. (|18fl . The only term in Vsb will be taken to be the simplest one given in Eq. 122fl . 
In this minimal picture, there is no symmetry breaking associated with the qqqq fields in M' . The symmetry breaking 
in the physical states (which contain two quark as well as four quark components) is due to the mixing terms which, 
as we have already seen in Eq. ifTfljl . can be invariant under SU(3)x SU(3)x \J(1) A . 

The method of treatment, as used earlier [3^| to discuss the model containing only the field M, is based on two 
generating equations which reflect the invariance of Vq under vector and axial vector transformations. Differentiating 
them once, relates two point vertices (masses) with one point vertices. Differentiating them twice relates three point 
vertices (trilinear couplings) with masses and so on. These are essentially tree level Ward identities. 

Under the infinitesimal vector and axial vector transformations we have, 

SvV - {Tr&W+^8 v S) + (<l>,S)->(<l/,S')} = ) 

6 A V = {Tr + ^aS) + (4>, S) - (0', S')} - -£„ (27) 
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wherein the non-zero value of the axial variation equation reflects the presence in Vo of the single U(1)a non-invariant 
term of Eq. (|18fl . Using Eqs. I|24[l and l|26|) as well as the arbitrariness of the variations Ey and Ea yields the matrix 
generating equations, 

{[0, ^] + [s, ^} + s) -> s')} = o, 

{[ ^< ^ ]+ - [s > w ]+ + ^ s) -> ^ s ' )} = 1[2Tr { *m - S W ] ~ 8C3Z ln (28) 

where, in addition, the form of Eq. (|18|) was used. To get constraints on the particle masses we will differentiate these 
equations once with respect to each of the four matrix fields: , S, S' and evaluate the equations in the ground 
state. Thus we also need the "minimum" condition, 

,dV dV SB ,_ ,dV ( BVsbs_ ( . (9q) 

Using our present choice of Eq. i|22|) as the only flavor symmetry breaker and Eq. l|12Jl , this becomes 

Q = 2A., (30) 

Now let us differentiate successively the vector generating equation with respect to S b and to S' a . This gives with 
the help of Eq. l|3U|) . the following two relations: 

d 2 V d 2 V 

(a a -a b ){ d2V \ ) + (&-&)( — , ) = 0. (31) 

The first of these equations relates the mass mixing transition with the unprimed scalar squared masses while the 
second of these relates the mass mixing transition with the primed scalar squared masses. It may be seen that 
information is obtained only for particles with different upper and lower SU(3) tensor indices. In the isospin invariant 
limit (where ot\ = 0.2 etc.), information will be obtained only for the kappa type particles (e.g. k + — Sf when mixing 
is neglected). If isospin violation information is inserted, information may be obtained also about the isovector scalars 
like a q (980) (which is represented by S 2 when mixing is neglected). Next, let us differentiate successively the axial 
vector generating equation with respect to 4> and to <fi'. It is neater to write the results first for the case when fields 
with different upper and lower tensor indices are involved: 



Next, let us write the corresponding equations for the case when the upper and lower tensor indices on each field 
are the same. 

d 2 V , Q , d 2 V , o , d 2 V , 8c 3 



MB 1 

d 2 V \ , a 1 d 2 V \- a ■ d 2 V 



'"a' )o t "<• ,.<"• I, g ""«"•'./ 



Note that the axial generating equation provides information on the masses of all the pseudoscalars. Further differ- 
entiations will relate a large number of trilinear and quadrilinear coupling constants to the meson masses and to the 
quark mass coefficients, A a . 

To fully characterize the system we will also require some knowledge of the axial vector and vector currents |3S| 
obtained by Noether's method: 

{j; xial ) b a = (a a + a b )d^ b a + (p a +p b )d^' b a + ---, 

(j; ect ° r )a = i{a a -a b )d ll Sl+i{(3 a -(3 b )d tl S ,b a + ---, (34) 
where the three dots stand for terms bilinear in the fields. 
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IV. PREDICTIONS FOR MASS SPECTRUM 



Here we consider the predictions for the mass spectrum of the model with the Lagrangian given in Eq. (|13|l . whose 
potential contains any SU(3)iX SU(3)jjX U(l)^ invariant terms whatsoever, amended with the SU(3)iX SU(3)jj 
but not U(1)a invariant term of Eq. (|18|l as well as the term, Eq.J22J) which transforms exactly like the QCD quark 
mass term. A characteristic feature is mixing between fields with the same quantum numbers. Specifically, there is 
information about mixing between tt and tt', between K and K' , between k and k' and among among the four 77 type 
(isosinglet) states. We will take these up in turn. Note that we will be working in the isotopic spin invariant limit 



A. The n — 7r' system 



For compactness let us denote, 



2Ax 

ai ' 
, d 2 V 



ai 



(35) 



Here we have introduced the total potential V = Vq + Vsb- However, since the second derivatives of Vsb vanish with 
our present choice of flavor symmetry breaker we may just use Vo- Substituting a = 1, b = 2 into both of Eqs. I|32(l 
enables us to write the (non-diagonal) matrix of squared 7r and tt' masses as: 



-\- z^y-jY 

^"kVt: 



Vtt 



(36) 



It is clear that is a measure of the mixing between 7r and tt' since the matrix becomes diagonal in the limit when 
z„ is set to zero. So we see that x^ would be the squared pion mass in the single M model and y v represents the 
squared mass of the "bare" tt' . Denoting the eigenvalues of this matrix by and m 2 ,, we read off the product and 
sum rules: 



mlml, = x^y„, 
+ y % (l + zl). 



(37) 



Assuming that the values of m* and m n i are known, the first of these equations expresses y^ in terms of x^ . Then 
the second of these equations also expresses z% in terms of x v . The value of x v is not known but its range is restricted 
to be, 



m„ < Xir < m„ 



This range may be derived by expressing z% in terms of x^ as mentioned and requiring z% > 0. 

The transformation between the diagonal fields (say ir + and ir' + ) and the original pion fields is defined as: 



" 7T+ " 




7T'+ 





COS! 

sin f 



— sine 
cos 8 



The explicit diagonalization gives an expression for the mixing angle 8„ : 



tan (20*) 



(38) 



(39) 



(40) 



which evidently is also known, up to a sign choice for z„, once x^ is specified. 

The mixing angle, 9^ can also be connected to the experimentally known value of the pion decay constant (i.e. the 
amplitude for the ir + meson to decay to two leptons). Substituting the expressions from Eq. [p?9"|) for 4>\ and (f>'\ in 
terms of the physical fields ir + and 7r' + into Eq. I|34|l yields, 



+ 



F v eV + 



[ax 
(ai 



02 J cos e/jr 
a 2 )sin8 7r - 



(A+^ 2 )sin0 W! 
(A + / 9 2 )cosc9 7r . 



(41) 
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We can then obtain a\ (in the isospin invariant limit) as, 



2(cos0 ff — ZirSinO-n) 

We then successively obtain A\ from the definition of x w , Eq. I|35|) and j3\ from the definition of z w , Eq, Q35[l. To sum 
up, specifying i x and the experimental quantities m^^m^i and determines all the other parameters of the 7r — n' 
system. 



B. The K - K' system 

The treatment of this system is almost exactly analogous to that of the tt — tt' system above when one defines the 
analogous variables, 

2{A 3 +A 1 ) 

x K = 



Vk 



a 3 + ai 

d 2 v x 



zk = ■ (43) 

Substituting a = 1, b = 3 into both of Eqs. (|32|l enables us to write the (non-diagonal) matrix of squared K and K' 
masses as: 



(M 2 )= XK + Z x yK ~~ ZK V K 
-zkVk Vk 



(44) 



This is observed to be identical to the expression for (M 2 ) in Eq. I|36|l when one simply substitutes everywhere K 
for 7r and K' for it' . Similarly, the four equations 1|37J) . I|38(l . Ij39(l and H4()|l continue to hold when one substitutes 
everywhere K for tt and K' for n' . Similarly, the K + decay constant, Fk is now defined from, 

(J^ al )l = F K d, t K + +F K ,d^K'+ + ■■■, 

F K = (ai +a 3 )cos9 K - + f3 3 )sm6 K , 

F K > = (a 1 +a 3 )sm9 K + (l3 1 +(33)cos9K.. (45) 
We can then obtain a 3 + a\ (in the isospin invariant limit) as, 

"3 + olx = — — — — — ■ (46) 

cos Ok — zkSiti Ok 

We then successively obtain A 3 + A\ from the definition of Xk and (3 3 + fix from the definition of zk- To sum up, 
specifying xk and the experimental quantities mK,mK' and Fk determines all the other parameters of the K — K' 
system. 



C. The k — k' system 

Again, we can treat this system in an exactly analogous way to the 7r — n' and K — K' cases if we define the 
analogous quantities: 

_ 2(^3-^1) 

d 2 v 

Vk K dS>\dS>\ h 

z, = AlA. (47) 
a 3 — ot\ 
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In this case, however, the vector generating equations in Eqs. (|31fl with the choices a = 1 and b — 3 are used. The 
transformation between the diagonal and original strange scalar fields is given by, 







k'+ 





cos V K — sin f 
sin 9 K cos 



si 



where the mixing angle is determined by the diagonalization: 



tan(29 K ) = 



-2y K z K 



(48) 



(49) 



We may define n "decay constants' as, 



F K = (a 3 - ai)cos6» K - (0 3 - 1 )sm6 K , 
F K > = (a 3 - ai)sin6>„, + (0 3 - /J^cos^, 



(50) 



although there is no direct experimental information available about them. 

Now let us consider the ir — n' , K — K 1 and k — k! systems together. Using the first two we can get all of 
A\, A 3 , a%, a 3 , 0%, 3 from the experimental masses of 7r, 7r', K, K' , the experimental decay constants F„, Fx and the 
assumed values of x„ and xk , as seen above. This means that x K and z K may be read off directly from Eqs. H47|) 
while y K can be found from the product rule m 2 K m 2 R , — x K y K if m K and m K t are furnished. Thus all the parameters of 
the k — k! system are known, given the input masses and the values of x^ and xk- However we have not yet made 
use of the sum rule analogous to the second of Eqs. 1(3*7)1. This provides another way to calculate z K so we get the 
consistency condition: 



1 N2 



x K {m 2 K 



Xk) 



(51) 



Since the quantities in this equation depend on both x^ and xk , the solution can determine the value of xk for each 



choice of x n . In other words, if x n is specified, the parameters of the tt 
all determined in the present model. 



the K-K' and the 



K ~ k systems are 



D. The r/ system 



This system is more complicated because, even in the isotopic spin invariant limit, there are four different 1 = 
pseudoscalars which can mix with each other. These may be put together as a column vector according to, 



$ 



1+<P 2 
I/O 



(52) 



The part of the Lagrangian describing the masses of the 1 = pseudoscalars is then: L = — (1/2)Qq (M 2 )Q , where 
(M 2 ) is a symmetric 4x4 matrix. Relations among the matrix elements follow by using both of Eqs. (|33l) . These 
connect the transition masses both to the "bare" unprimed particle masses and to the "bare" primed particle masses. 
The use of isospin invariance relations like the ones given in Appendix B may also be useful. Eventually, the matrix 
elements of (M 2 ) depend on four new quantities in addition to the ones appearing in the above three subsystems. 
The resulting matrix elements are listed below: 



K 2 )n 



2Ai 
+2 



16c 3 f3lm 2 m 2 n , 



2A x ai 



(- 



d 2 V 



aia 3 



\/2Ai 



Ot3 



8V2c 3 Plmlml, ( 2y/20( 



d 2 v 



l\2< 



\ ai"3 J (dcjj'i) 



d 2 V , 

1\2< 



d 2 V 



d 2 V 
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(M 2 ) 



K 2 ) 

( M v) 



2\ 

14 



2^ 

22 



(M 2 ) 



K 2 ) 

K) 



21 




r) 2 V 

(M 2 ) A , = ( - „ ) (53) 

The four new quantities are C3, discussed earlier, and the "bare" primed squared masses: 

, d 2 V d 2 V d 2 v 

Wl) 2 '' Wi^l Wi^T ' 

These four quantities may be found by inputing the masses of four isosinglet pseudoscalars. The net result is that all 
four systems discussed will be completely described if all the experimental masses and the decay constants, F^Fx 
are specified together with an assumed value for x„. 



V. COMPARISON WITH EXPERIMENT AND DISCUSSION 



In the preceding section we gave the tree level formulas resulting from the M — M' model with any SU(3)lX 
SU(3)^ x U(l)^4 invariant terms together with a single "instanton" type term which mocks up the U(l)^ anomaly and 
the simplest structure which mocks up the quark mass terms. Isotopic spin invariance was also assumed. Information 
is provided for only the pseudoscalar nonets and the strange scalar particles. Information about the scalar isotriplets 
can be obtained by including isospin violation effects while information about the scalar isosinglets requires either 
assuming some specific form for the invariant interaction terms or computing other physical quantities. These will be 
discussed elsewhere. Now we will input the experimental masses to try to learn what the model has to say about the 
quark structure of the various mesons being described. In particular we are interested in the mixing angles like (9^, 
governing admixtures of qq and qqqq in the physical states and the four quark "condensate" strengths, f3 a which are 
associated with this mixing in the present model. 

The well known lowest pseudoscalar nonet masses and decay constants will be taken, for definiteness (considering 
the ambiguity as to which member of a non trivial isospin multiplet to choose), to be: 

m,, = 0.137 GeV , m K = 0.496 GeV, 
m v = 0.548 GeV , rry = 0.958 GeV, 

F n = 0.131 GeV , F K = 0.160 GeV. (55) 

Next, let us consider what are the suitable experimental inputs for the masses of the excited mesons, w',K'.k' 
and for the k meson itself. In the latest Review of particle properties 01 there are two dotted (i.e. considered 
established) candidates for excited pions below 2 GeV: the 7r(1300) and the 7r(1800). These particles could have 
four quark components and/or radially excited two quark components. In fact, judging from an investigation of 
excited baryons |48|. it is likely that both types are present. Clearly, however, for our present investigation it seems 
reasonable to assume that the four quark component is the dominant one and to choose the lower mass object as the 
more suitable one. Similarly there are two undotted (non established) excited kaon candidates: the K(1460) and the 
K(1830). We will again choose the lower value. As candidates for an excited strange scalar there is a dotted i"CQ(1430) 
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and an undotted I^q (1950) and we again choose the lower value. In the case of the low mass strange scalar there is an 
undotted if q (800) candidate, which we will interpret, with the help of 15], to be closer to 900 MeV. We summarize 
these choices: 

rrv = 1-30 GeV , m K > = 1.46 GeV, 

m K = 0.90 GeV , m K > = 1.42 GeV. (56) 

For the excited r\ type pseudoscalar particles the Review of particle properties lists, below 2 GeV, the possible 
masses (all in GeV): 

1.294, 1.410, 1.476, 1.760. (57) 

The first three of these are dotted but the fourth is undotted. Here it seems more difficult to a priori choose which 
are most relevant so we shall study all possible pairings in a systematic way. 

First let us discuss the tt — n', K — K' and k— k' systems. After using the inputs of Eq. (|55J) . all features of these 
systems in our model will, as already discussed, be determined by specifying x n . Table[I]shows the predicted physical 
parameters for three values of x w . For orientation we note that in the chiral model with a single field, M one has 

a x -> F«/2 = 0.0655 GeV , a 3 -> F K - a x = 0.0945 GeV, 

A x -► ymj = 6.15 x 10~ 4 GeV 3 , A 3 -> -ym^ = 0.01866 GeV 3 , 

fa -> , fa - 0. (58) 

The single M model corresponds to the choice x v = mj. Increasing x„ has the effect of increasing the admixture 
of the "four quark" field component in the physical pion. The "quark mass ratio", A3 /Ax = 30.3 in the single M 
model is not very different from the value of 31.2 obtained using the values in the i T = 0.019 Gev 2 column. The qq 
meson condensates ax and a 3 are also very similar. Of course the "four quark" meson condensates fix amd fa are zero 
without M' . Despite the similarities, the 6.4° mixing angle already corresponds to about an 11 percent "four quark" 
admixture in the physical pion wave function. Considering that the accuracy of current algebra predictions for low 
energy pion physics is roughly ten percent, it seems that this choice of x„ is the most plausible one. One sees from 
the second and third columns that relatively small increases in x n lead to large increases in four quark admixture for 
the pion and the kaon. Interestingly, the behavior of the four quark admixture in the strange scalar meson k is quite 
different. When the pseudoscalars are closer to pure "two quark" states in the model the scalar has a large four quark 
admixture (34.1°, with the choice of x„ in the first column). Thus the result is consistent with having a fairly large 
four quark component in the light scalars. 





3^=0.019 (GeV 2 ) 


av=0.021 (GeV 2 ) 


2^=0.022 (GeV 2 ) 


On (deg.) 


-6.4 


-19.1 


-22.7 


0k (deg.) 


-11.2 


-22.9 


-26.2 


On (deg.) 


34.1 


28.1 


26.5 


Ai(GeV 3 ) 


6.19 x 1(T 4 


6.51 x 1(T 4 


6.66 x 10" 4 


A 3 (GeV 3 ) 


1.94 x 10" 2 


2.07 x 10~ 2 


2.12 x 10~ 2 


ai (GeV) 


6.51 x 1(T 2 


6.20 x 10~ 2 


6.06 x 10~ 2 


a 3 (GeV) 


9.24 x 10" 2 


8.83 x 10~ 2 


8.69 x 10 -2 


Pi (GeV) 


7.18 x 1(T 3 


2.12 x 10~ 2 


2.50 x 10~ 2 


03 (GeV) 


2.03 x 10~ 2 


3.38 x 10" 2 


3.74 x 10" 2 



TABLE I: 9k, 0k and 6 K are respectively the "four quark" admixtures in the n,K and k states. Ai,As represent the quark 
mass parameters while ax, 03 and Px,Pz represent respectively the two and four quark condensate strengths. These are plotted 
as functions of the assumed "bare" pion squared mass, Xn- 

The analogs of the two quark condensates ax = Q?2 and 0:3 are approximately equal, in agreement with the usual 
assumption that the vacuum is approximately SU(3) symmetric. The analogs of the four quark condensates in this 
model are roughly an order of magnitude smaller than the similarly normalized two quark condensates. They are 
furthermore seen to deviate appreciably from SU(3) symmetry. It should be noted, as discussed in ref, |16| for example, 
that the tensor indices for the primed mesons really correspond to "dual quark" or diquark indices in accordance with, 



Qa ~ tab c q b q c - 



(59) 
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Thus in terms of the usual quarks, 

Pi r~j (ddss), 02 ~ (uuss), fli ~ (dduu). (60) 

Now consider the mixing of the four r\ type fields in the model. The basis is given in Ea. (|52|) while the elements 
of the 4x4 mass squared matrix are given in Eq. i|53[l . The orthogonal transformation matrix, K which relates the 
mass eigenstate fields, $ to the original ones is defined by 

$ = K§. (61) 

As discussed in the previous section, there are, after using the symmetry information, four new unknown parameters 
characterizing the r\ system. Thus taking the four mass eigenvalues from experiment could in principle determine, 
together with results from the 7i — %' , K — K' and k — k' systems, everything about the r\ system for a given value 
of x n . However there is no guarantee that there will be an exact solution for all choices of experimental parameters. 
This is the case, in fact, so we will search numerically for a choice of "theoretical" masses which will best fit the 
experimental inputs. The criterion for goodness of fit will be taken to be the smallness of the quantity: 

x ^Y,\ m T v --^ co M^T v -- (62) 

i 

As shown in Eq. I|57|) . there are three established candidates and one not yet established candidate below 2 GeV for 
the two excited rj states. This yields six possible scenarios for choosing them. The quantity x f° r each choice is shown 
in Table ITU for three values of the parameter x^. It may be observed that the fits typically get worse with increasing 
x^, so it is reasonable to consider the choice 0.019 GeV 2 for this quantity as we did previously. The smallest values of 
X are found for scenarios 5 and 6. However these both involve the 77(1760) state which is the one not yet established. 
The smallest value of x using only established states is scenario 2. This case corresponds to an exact fit with eta type 
masses in GeV (experimental values in parentheses for comparison): 

0.533(0.548), 0.963(0.958), 

1.327(1.294), 1.716(1.476). (63) 



Scenario 


2^=0.019 (GeV 2 ) 


xv=0.021 (GeV 2 ) 


3^=0.022 (GeV 2 ) 


1 


{t/(1295), 7/(1405)} 


6.23 x 10~ 2 


3.99 x 10" 1 


5.08 x 10" 1 


2 


{77(1295)^(1475)} 


2.85 x 10~ 2 


3.39 x 10" 1 


4.44 x 10" 1 


3 


{77(1295), 77(1760)} 


2.35 x 10 -2 


1.37 x 10 _1 


2.28 x 10 _1 


4 


{7/(1405), 77(1475)} 


8.28 x 10~ 2 


3.63 x 10" 1 


4.49 x 10 _1 


5 


{7/(1405), 7/(1760)} 


1.50 x 10~ 2 


1.62 x 10 _1 


2.38 x 10 _1 


6 


{7/(1475), 7/(1760)} 


2.84 x 10~ 2 


1.78 x 10 _1 


2.68 x 10 _1 



TABLE II: A goodness of fit quantity, % = ]Ci l m i XP ~~ m* heo '| /m" xp ', where the TTii are the four mass eigenvalues of the rj 
type fields, is given for 6 possible scenarios and for three values of x v . Each scenario corresponds to a choice of 77 type fields 
including the 7/(548) and the 7/(958) as well as the two listed in the left hand column. 



The detailed content of all the 77 mass eigenstates can be read off from the matrix K 1 . For scenario 2 we have, 



K~ 



-0.570 0.750 -0.023 0.333 

-0.329 -0.573 0.142 0.737 

0.704 0.267 -0.309 0.581 

0.267 0.192 0.940 0. 



(64) 



Thus, in the present model there is an 89 percent probability ((K _1 )f 1 + {K^ 1 }^) that the t/(548) is a quark-antiquark 
state and an eleven percent probability that it is a four quark state. As expected, the 77(548) is most likely to be in 
an ss state. In the case of the 77(958), there is a 44 percent probability for it to be in a quark antiquark state. There 
is a 54 percent probability for it to be in the four quark state (f>' 3 . This situation has some plausibility since in terms 
of ordinary quarks, the latter state has the content uudd and it should be most energetically favorable to bind a four 
quark state made without strange quarks. 



13 



The other scenarios which don't employ the unconfirmed ?y(1760) state (numbers 1 and 4) have contents very similar 
to the one in Eq. (|fc>4|) . On the other hand the three scenarios employing the ^(1760) have a rather different content, 
which seems unusual: scenarios 3, 5 and 6 make the 77(958) almost completely <f>' 3 . 

In scenario 2, which seems the most reasonable choice, we notice that the ry(1295) has a 43 percent probability of 
being in a four quark state while the 77(1475) has an 89 percent probability of being in a four quark state. To sum up, 
the value x„ =0.019 GeV 2 leads to fairly small four quark content in the light pseudoscalars- tt, K, r\ at the same time 
that the light scalar k has an appreciable four quark component. The "excited" ?7's are predominantly four quark 
states. The 77(960) is mainly two quark in content but has a non trivial four quark piece. 

The results obtained here provide supporting evidence for the feature, illustrated in the first treatment of this model 
|29) . that the lightest scalars, unlike the lightest pseudoscalars, have appreciable four quark components. That model 
neglected quark masses and used the simplified choice of terms shown in Eq. I|14|) . The more recent treatment of 
ref. includes two additional invariant terms beyond those in Eq. 1)14(1 (although not all the renormalizable terms 
shown in Appendix A) as well as four types of quark mass splitting terms. Our results for the present treatment, 
where quark masses are included and which holds for any possible SU(3)lX SU(3)_rX U{1)a conserving terms, are 
also in qualitative agreement for the ?r-type, K-type, 77-type and K-type states with that treatment. Roughly, this may 
be expected since the present approach includes any choice of invariant terms. However, we only used here the single 
quark mass splitting term of Eq. ((22(1 . Thus the results seem qualitatively robust with respect to the treatment of 
the mass splittings. 

An interesting feature of our model is the presence of "four quark" condensates as signaled by the non-zero values 
of the p a - To make a rough estimate of what this corresponds to in quark language we proceed as follows. In ref. fltl ] 
it was pointed out that the mass formulas of the single M linear sigma model could be transformed to the "current 
algebra" ones 49] by the replacements: 

A a = m a A 2 , a* = (65) 

where the m a are the ("current" type) quark masses and A is the QCD scale factor. Taking A\ = 6.19 x 10~ 4 GeV 2 
from the left column of Table HI and mi « 5 MeV we get A w 0.35 GeV (and (q a q a ) ~ -0.016 GeV 3 ). In the case 
of the four quark condensate, as one sees from the discussion in the Introduction, there are several ways to couple 
the four quarks together to make scalars. We are assuming that one such way has been selected. For that case, it is 
reasonable to expect, on dimensional grounds, that 

\(ddss)\ ~ A 5 /3i «4x 10" 5 GeV 6 . (66) 

In comparing the scalar masses with experiment there are expected to be, as discussed in the first four sections of ref. 
[2^|, non- negligible corrections due to the use of unitary models for the pseudoscalar- pseudoscalar scattering based 
on this Lagrangian. We plan to report on this elsewhere. This should also enable us to study the isosinglet scalar 
masses. For both isosinglet scalars and pseudoscalars, the inclusion of possible glueball states is another interesting 
topic we plan to pursue. The additional symmetry breaking terms like those in Eqs. ((19(1 and 1(23(1 seem also to be 
worth investigating. 
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APPENDIX A: RENORMALIZABLE MODEL 

The twenty one SU(3)x SU(3) invariant renormalizable terms in Vq of Eq. 1(1311 which can be made out of M and 
M' are: 

V = - C2Tr(MM t ) + £ 3 (detM + h.c.)+c^Tr(MA7 rt MM t ) + c^ (Tv(MM r )) 2 

+ d 2 TrfM'M'tjTrfMM't) + d 3 (detAf + h.c.) + d\ Tt{M'M'^M'M^) + d\ (Ti(M' M'^)) 2 
+ e 2 (Tr(MM't) +h.c.) 

+ e a 3 {e ahc e de t ' M a d M h e Mf + h.c.) + e\{e abc e de f M%M*Mf + h.c.) 
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+ e\ Tr(MM f M'M n ) + e\ Tr(MM'^M'M^) 

+ e\ [Tr(MM'*MM't) + h.c] + ej [Tr(MM^ MM'^) + h.c.] + e\ [Tr(M'M' t M'M t ) + h.c] 

+ e{ Tr(MM t )Tr(M'M' t ) + el Tr(MM' t )Tr(M'M t ) + e£ [(Ti^M'M't)) 2 + h.c] 

+ 4 [Tr(MM t )Tr(MM' t ) + h.c] + [Tr(Af'M' t )Tr(M'M t ) + h.c]. (Al) 

Notice that among these terms, those with the coefficients C2, c% 1 c\, d%, d\, e 3 , ef , e^, e{, ef and are U(V)a invariant. 
It also may be of some interest to write down the twenty one renormalizable terms, linear in the matrix A, which 
transform like the QCD quark mass terms under SU(3)x SU(3). Again, for this listing, the U(l)^ transformation 
property of the mass terms in the fundamental QCD Lagrangian is respected only for the terms with the coefficients 

ki, &3, &4, kg, fcll,fcl2, kl7, k 2 \. 



V SB = + ki[Tr(AM) + h.c.]+k 2 [Tr(AM') + h.c] 

+ k 3 [Tr(AMM^M) + h.c] + fc 4 [Tr( J 4MAf' t A/ / ) + h.c] 

+ k 5 [Tr(AM M^M') + h.c] + k 6 [Tr(AMM n M) + h.c] 

+ k 7 [Tr(AM'M'^M') + h.c] + ks [Tr(AM'M^M) + h.c] 

+ kg [Tr(AM'M^M) + h.c] + k w [Ti(AM' 'Aft AT') + h.c] 

+ fcu[Tr(AM) +h.c.]Tr(A/A/ t ) 

+ k 12 [Tr(AM) + h.c]Tr(A/'Af' t ) 

+ k 13 [Tr(AM)Tr(MM^) + h.c] + fci 4 [Tr(AAf)Tr(Af 'M f ) + h.c] 

+ k 15 [Tr(AM') +h.c]Tr(AfAT t ) 

+ k 16 [Tr(AM') +h.c.]Tr(Af'M' t ) 

+ k 17 [Tr(AM')Tr(MM' r ) + h.c] + k 18 [Tr(AM')Tr(M' M r ) + h.c] 

+ k 19 A b a e bcd e ae fM*Mf + h.c. 

+ k 20 A b a e bcd e ae f M' c e M' d f + h.c. 

+ k 21 A b a e bcd e ae fM^M' d f +h.c. (A2) 



APPENDIX B: SOME ISOSPIN RELATIONS 



We give examples of relations which follow from isotopic spin invariance: 



d 2 v x 


, d 2 v 




'd<t>\d<f>\ 


d 2 V x 


, d 2 v 


d<t> 2 2 d<j>l } ~ 






1 d(f>{d(f> 3 3 


d 2 V 


, d 2 v 




'd<t>\d<t>\ 



d 2 V 

(Bl) 



"d<l>{d<l>i 



Similar relations hold when V is differentiated with respect to two primed fields and with respect to one primed and 
one unprimed field. 
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